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ABSTRACT 

We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic 
spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located 
at distances below 5.5 kpc, identified from the ATLASGAL 870yum data. We define three evolutionary classes within this sample: 
starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with Hii regions. We find that the 
NiHi) probability density functions (A-PDFs) of these three classes of objects are clearly different: the A-PDFs of starless clumps are 
narrowest and close to log-normal in shape, while star-forming clouds and H ii regions exhibit a power-law shape over a wide range 
of column densities and log-normal-like components only at low column densities. We use the A-PDFs to estimate the evolutionary 
time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the 
A-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds 
contain greater relative amounts of dense gas across all evolutionary classes. 
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1. Introduction 


Molecular clouds (MCs) are the densest regions of the interstel¬ 
lar medium and the birth sites of stars. Nevertheless, despite this 
important role in star formation, key aspects of MC evolution 
remain unclear: What are the key parameters in determining the 
star-forming activity of MCs? How do these parameters change 
with MC evolution? The column density distribution of MCs 
has been found to be sensitive t o the relevant physical processes 
(iHennebelle & Falgarone|[20T^ . The study of the density struc¬ 
ture of clouds that are at different evolutionary stages can there¬ 
fore help to understand which physical processes are dominating 
the cloud structure at those stages. 

Column density probability density functions (A-PDFs) are 
useful tools for inferring the role of different physical pro¬ 
cesses in shaping the structure of molecular clouds. Observa¬ 
tions have shown that non-star-forming molecular clouds show 
hottom-heavj^N-V\yFs, w hile the star-forming molecular clouds 
show top-heav^N-V\y¥s (iKainulainen et al.ll2009[l2011~b[l20l4 
Kainulainen & Tan|[2013l: ISchneider et al.ll2013h . It is generally 
accepted that the top-heavy A-PDFs are well described by a 
power-law function in their high-column density regimes. The 
description of the shapes of the low-column density regimes of 
both kinds of A-PDFs is still a matter of debate. The papers cited 
above describe the lo w-column de n sity re gime s as log-normal 
functio ns. In contrast, lAlves et ^ (l2Q14l) and iLombardi et al.l 
(l2Q15l) argue that a power-law function fits the observed A-PDFs 
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^ Most of their mass is in low-column density material. 

^ They have a significant amount of mass enclosed in high-column 
density regions. 


throughout their range. The origin of these differences is cur¬ 
rently unclear. 

Simulations predict that turbulence-dominated gas devel¬ 
ops a log-normal A-PDF (iFederrath & Klesse^ |2013|) : such 
a form is predicted for the volume density PDF (hereafter 
p-PDF) of isot hermal, supersonic turbulent, and non-self- 
gravitating gas (IVazQuez-Semadenil 1199^ iPadoan et al.l \l9^ 



Scalp et al. H998I: lOstriker et al.l I 2 OOII: JPadoan & Nordluni 




iBallesteros-Paredes et al.l l201 ll: IFederrath & Kless^ 
Log-normal p-PDFs can, however, result also from 
processes other than supersonic turbulence such as gravity op¬ 
posed only by ther mal-pressure forces or gravitationally-driven 
ambipolar diffusion (iTassis et al.ll201Ql) . 

The log-normal A-PDF is defined as: 


p(v;p,(T^) = 


1 


: 


exp 


-(s-pY 

2(T? 


( 1 ) 


where 5- = In (Ay/Ay) is the mean-normalized visual extinc¬ 
tion (tracer of column density, see Section m, and p and 
CTs are respectively the mean and standard deviation of the 
distribution. The log-normal component that is used to de¬ 
scribe low column densities has t ypically the width of = 
0.3 - 0.4 (IKainulainen et al.l[20091) . It has been suggested that 
the determination of the width can be affected by issues such 
as un related dust emission along the line of sight to the 
cloud (ISchneider et al.ll20T^ . Practically all star-forming clouds 
in the Solar neighborhood show an excess to this component at 
higher column de nsities, following a power - law, or a wider log- 
nornial function ([Kainulainen et al.l [20091: IKainulainen & TanI 
l2013l:[^hneider et al.ll2(313l). especially reflecting their ongoing 

star formation activity (IKainulainen et al .112014 ISadavov et al.l 
l20l4 IStutz & Kainulainenll2015l) . Such behavior is suggested 
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by the predictions t hat develope top-heavy p-PDFs for self- 
gravitating systenis (iKlessenI lioort Iklessen & BurkertI I200QI: 


gravitating systenis (IKlessenI IZUUU IKle: 

iKritsuk et al.ff2Ql it iGirichidis et al.|[20T4l) 

Another interesting measure of the density structure of 
molecular clouds is the dense gas mass fraction (DGMF) that 
describes the mass enclosed by regions with M{Av > Ay), rela¬ 
tive to the total mass of the cloud, Mtot- 


dM' = 


M(Av > A'y) 
Mtot 


( 2 ) 


The DGMF has been recently linked t o the s tar-fo rming rates 
of molecula r clouds: iHeiderman et al.l (l201Ql) and iLada et aP 
(|201QL[20T^ showed, using samples of nearby molecular clouds 
and external galaxies, that there is a relation between the mean 
star-forming rate (SFR) surface density (Ssfr) and the mean 
mass surface density (Smass) of MCs: Esfr oc foc'^mass, where 
/dg = M{Av:^.0mag) ^ ^ Sample of eight molecu¬ 

lar clouds within 1 kpc, a c orrelation Esfr ^ ^mass reported 
by iGutermuth et al.l (l201ll) . Combining these two results sug¬ 
gests /dg ^ ^mass- 

Despite their utility, a complete, global understanding of 
the A-PDFs and DGMFs of molecular clouds is still missing. 
One of the main problems arises from the fact that the dynamic 
ranges of different observational techniques sample the A-PDFs 
and DGMFs differently. Previous works have employed various 
methods: CO line emission only samples A-PDFs between Ay 
3-8 mag (iGoodman et al .11200% . NIR extinction traces column 
density at wider, but still narrow , dynamic range. Ay ^ 1 - 
25 mag (iLombardi & Alvesll200ll) . iKainulainen & TanI (l2Q13l) 
and iKainulainen et al.l ( 2013al) used a novel extinction technique 
that combines NIR and MIR data, considerabl y increasing the 
obsery able dynamic range. Ay = 3 - 100 mag. [Schneider"^ III 
(1201 3l) and iLombardi et al.l (1201 5l) used Herschel FIR data to 
sample A-PDFs at Ay < 100 mag. 

Another observational hindrance arises from the limited spa¬ 
tial resolution of observations. The high-column densities typi¬ 
cally correspond to small spatial scales in molecular clouds; to 
probe the A-PDFs at high-column densities requires spatial res¬ 
olution that approaches the scale of dense cores in the clouds 
(~0.1 pc). Herschel reaches resolution of ~36" that corresponds 
to 0.17 pc at 1 kpc distance. Extinction mapping using both NIR 
and MIR wavelengths can reach arcsecond-scale resolution, but 
only a bout ten clouds have been studied so far with that tech¬ 
nique (IKainulainen & Tanll2013l:[^inulainen et al.ll2013E) . 

Born out of the observational limitations above, the most im¬ 
portant weakness in previous studies is that they only analyze 
relatively nearby molecular clouds (J < 1.5 kpc). Thus, they 
probe only a very limited range of Galatic environments, which 
prohibits the development of a global picture of the factors that 
control A-PDFs across different Galactic environments. Extend¬ 
ing A-PDF studies to larger distances is imperative for three 
principal reasons. First, studying the more massive and distant 
MCs will allow us to sample the entire MC mass range present 
in the Galaxy. Second, larger numbers of MCs over all masses 
provide statistically meaningful samples. Finally, extending to 
larger distances is necessary to study the possible effect of the 

Galactic structure on the mass distribution statistics^^ _ 

_ In this paper, we emplo y the ATLAS GAL (ISchuller et al.l 

120091: ICsengeri et all 12(^1) ' survey to study a large sample 
of molecular clouds in the Galaxy. The ATLAS GAL survey 
traces submillimeter dust emission at 870yum. Submillimeter 
dust emission is an optically thin tracer of interstellar dust, and 
hence a direct tracer of gas if a canonical dust-gas mass ratio 


is assumed. The submillimeter observing technique employed in 
the ATLASGAL survey filters out diffuse emission on spatial 
scales greater than 2.5', hence making the survey most sensi¬ 
tive to the densest material of the interstellar medium in which 
star formation occurs. With this data set we can observe the cold 
dense interiors of molecular clouds in both the near and far sides 
of the Galactic plane. With an angular resolution of 19.2", AT¬ 
LASGAL improves by almost a factor of two the resolution of 
Herschel observations, thus providing a more detailed view of 
the dense material inside molecular clouds. We will use this data 
sample to study the A-PDFs and DGMFs of molecular clouds at 
different evolutionary classes. 


2. Data and methods 

We used continuum maps at 870yum from ATLASGAL to iden¬ 
tify MCs in the Galactic plane region between I e [9°, 21°] and 
\h\ < 1°, where the rms of the survey is 50mJy/beam. We se¬ 
lected this area, because extensive auxiliary data sets were avail¬ 
able for it, a nd specifically, st a rless clumps have already been 
identified bv iTackenberg et all (l2012h . We classified the identi¬ 
fied molecular cloud regions in three groups based on their evo¬ 
lutionary classes: starless clumps (SLCs), star-forming clouds 
(SFCs), and H ii regions. In the following, we describe how each 
class was defined and how we estimated the distance to each re¬ 
gion. 


2.1. Source selection 


We identified molecular cloud regions based primarily on AT¬ 
LASGAL dust emission data. As a first step, we defined objects 
from ATLASGAL data simply by using 3cr emission contours 
(0.15 Jy/beam) to define the region boundaries. Then, we used 
distances available in literature (see Sect. 12.21) to group together 
neighbouring objects located at similar distances (within the as¬ 
sumed distance uncertainty of 0.5 kpc), i.e., those that are likely 
associated with the same molecular cloud. As a next step, we 
expanded the boundaries of the regions down to their Icr level 
in the cases in which they show close contours at Icr. Finally, 
each region created in this manner was classified either as a 
SLC, SEC, or H ii region using information about their stellar 
content available in literature. An example of the region defini¬ 
tion is shown in Fig.[T](see also Appendix[0). We identify a total 
of 615 regions, 330 of them with known distances and classified 
either as SLC, SEC, or Hii regions (Fig.O. Throughout this pa¬ 
per we refer to each of the ellipses shown in Fig. [T] with the term 
region. In the following we explain the definition of the three 
evolutionary classes in detail. 

HII regions are defined as regions ho sting previously cata¬ 
loged H ii regions . We us e d the catalog u es IWood & Churchwef 


(11989 ). iLockman (Il989l). iGarav et ^ (Il993l). iBronfman et aL. 
( 1996). iLockman et all (Il996l) . iForster & Caswelll (l200Ql) . 


and lUrquhart et al.l (l2013l) . We identified 114 H ii regions in the 


considered area. Distances are known for 84 of them (74%). Two 
thirds (57) of the H ii regions with distance estimates lie at near 
distances (J < 5.5 kpc). If we assume the same distribution for 
the 30 HII regions with unknown distances, 20 of them would 
be located at near distances. Nevertheless, we exclude these re¬ 
gions from our analysis. We summarize the number of regions 
with and without known distances in Table [T] (see Sect. 12.21) . 

The star-forming clouds (SFCs) are defined as the subset 
of regions devoid of H ii regions but containing young stel¬ 
lar objects (YSOs) and protostars. Here the presence of YSOs 
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and protostars is assumed to be a clear indication of ongo¬ 
ing star fo rmation. For this pur pose, we used t h e YSO cat¬ 
alogues of iDunham et al.l (l2Qlll) and iToth et all (l2Q14l) . The 
former search signs of active star formatio n in the Bolocam 
Galactic Plane Survey (I Aguirre et al. 20 lH BGPS) usin g the 
GLIMPSE Red Source catalogue ([Robitaille et al.l l200"^ . the 
EGO catalogue (ICvganowski et al.l l2008l) . and the RMS cata¬ 
logue (iLumsden et al.ll2013l) . They found 1341 YSOs in the area 
I € r9 °,2l°1 and IZ^I ^ 0.5° a nd it is >98% com plete at the 0.4 Jy 
level (iDunham et al.ll201 ih . iToth et all (l20l4 present a catalog 
of 44001 YSO candidates, 2138 in the area I e [9°,21°] and 
\b\ < 1°, with a reliability of 90% in the YSO classification. All 
the regions showing spatially coincident YSOs were classified as 
SECs. We only require one YSO to classify a region as SEC, but 
our SECs have more than one. The probability of classifying a 
SEC as a region without YSOs due to completeness issues in the 
YSOs catalogues is therefore very low. We identified 184 SECs, 
126 of them with known distances. The 80% (99) of the SECs 
with known distances lie at J < 5.5 kpc and are therefore stud¬ 
ied in this paper. Assuming the same SEC distribution for the 
SECs with unknown distances, we estimate that the 80% (46) of 
the SECs with no distance estimates would be located at near 
distances. 


E inally, we adopted the starless clump catalog 
from iTackenberg et al.l (l2012h to define our sample of SLCs. 
They present a SLC sample with peak column densities 
N > 10^^ cm“^. The properties of this SLC sample were 
specifically c hosen in order to detect potential high-mass star 
progenitors. ITackenberg et all (l2012h used uniform criteria to 
classify their SLC san iple: absence of GLIMPS E and/or 24 yum 
MIPSGAL sources. ITackenberg et al.l (l2012l) identified 120 
SLCs with known distance^ in the Galactic plane area studied. 
All SLCs are located inside our previously defined H ii regions 
or SECs (see Eig.[T]). 

We note a caveat in the above evolutionary class definition 
scheme. Our scheme makes an effort to capture the dominant 
evolutionary phase of the region, but it is clear that not all the 
regions are straightforward to classify. In principle, the distinc¬ 
tion between H ii regions and SECs is well defined; it depends on 
whether the regions host an H ii region or not. However, eight re¬ 
gions harbor only UCH ii regions whose extent is tiny compared 
to the full extent of those regions (#34, #54, #192, #195, #233, 
#246, #247 and #390). Since our aim is to capture the dominant 
evolutionary phase, we classified these regions as SECs. 

We also note that our evolutionary class definition is based 
only on the stellar content of the regions. The SLCs exhibit 
no indications of star-forming activity, SECs have star-forming 
sources, H ii region have formed massive stars. However, we 
emphasize that we cannot assume that all the SLCs will defi¬ 
nitely form stars. Similarly, we cannot assume that all the star¬ 
forming content within SECs will become massive enough to 
create H ii regions, although some of them will. Therefore we 
do not aim to draw a sequential evolutionary link between these 
three classes of regions. Instead, the estimated time-scales for 
each class instead aim to identify independent evolutionary time- 
scales for each observational class. 


^ We adopt only regions with solved kinematic distance ambiguity 
(KDA) as sources with known distances. 


2.2. Distance estimates and convolution to a common spatial 
resolution 


We adopted distances to each regio n from literature. The 
two m ain lit erature source s used were lEllsworth-Bowers et al.l 
(12013l) and IWienen et al.l (l2Q12h . The former catalog mea¬ 
sures kinematic distances of molecular clumps identified 
with sub-mm dust emission. They solve the kinematic dis¬ 
tance ambiguity (KDA) using Bayesian distance probability 
density functions. They use previous data sets to establish 
the prior distance probabilities to be used in the Bayesian 
analysis. This method has a 92% agreement with Galac¬ 
tic Ring Survey based dista nces. In total, 68 out of 33 0 
regions have counte rparts in lEllsworth-Bowers et al.l (l2013l) . 
IWienen et al.l (l2012h measured the kinematic distances to 
dense clumps in the ATLASGAL survey using ammonia ob¬ 
servations. We obtained distance estimates for 80 regions 
from this catalog. We also used other catalogs ba s ed on 


kinematic distances (IWalsh et al.l 1 199^ l^iinon et al.l 1200^ 


Rudolph et al.l 120061: [Roman-Duval et al.l 120091: lUrquhart et al.l 


20 131: ITackenberg et al .1120121). and in alhree-di mensional model 


of interstellar extinction (iMarshall et al.ll2009l) . A detailed dis¬ 
cussion on the methods for distance estimates is beyond the 
scope of this paper. We therefore refer to the cited papers for 
a detailed discussion on them. Table [2] shows the number of dis¬ 
tance estimates adopted from each literature source. In regions 
with more than one distance estimate, we estimated the distance 
averaging the different values. Eor all but six of the studied re¬ 
gions (~ 96%) the distance ambiguity was solved in at least 
one of the cited papers. Regions with different KDA solutions 
in literature (i.e. with several clouds along the same line-of- 
sight) were removed from our sample to avoid line-of-sight con¬ 
tamination. Eor the remaining six regions we used maps from 
the GLIMPSE and MIPS surveys to search for dark shadows 
against background emission (e.g. Stutz et al. 2009; Ragan et 
al. 2012). The near distance was adopted for regions associated 
with IRDCs. 


Since all the SLCs of our sample are embedded in 
Hii regions or SECs (see Section |2T]), we compared the dis¬ 
tance estimates for the SLCs and for their hosting regions. In 
every but one case, the distance estimates of the SLCs and their 
hosting SECs or H ii regions were in good agreement. In the 
onl y inconsistent case, the SLC was located at the far distance 
in (ITackenberg et al.ll2012l) and its hosting SEC was located at 
the near distance. Since the KDA solutions of the SLC and its 
hosting SEC differ, we removed out both regions from the final 
sample (see also previous paragraph). 

Eigure |2] shows the distance distribution of our sample. A 
vast majority (~ 80%) of our regions is located within 5 kpc dis¬ 
tance. There is a gap between 6 and 10 kpc, coinciding with the 
centra l hole of the Galactic molecular ring (ISolomon & Rivolol 
ll989lPI . At the far side of the Galaxy, there are three density 
enhancements that coincide with the Sagittarius, Norma, and 
Perseus spiral arms (Pig. ID). 

We study only regions within 5 kpc since the highest source 
density of our sample is located there. Assuming an error in 
distance determination of about 0.5 kpc, we also included re¬ 
gions located between 5 kpc and 5.5 kpc. We convolved the 
ATLASGAL data of all closer regions to a common 5 kpc 
distance resolution using a Gaussian kernel of FWHM = 


^ We note t hat the existence of the G alactic Ring has recently been 
questioned bv iDobbs & Burkerl (120121) . who proposed a two symmetric 
spiral arm pattern for the Milky Way as an explanation of observations. 
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Fig. 1. MIPSGAL 24 pm map of the Galactic plane between 9 deg < / < 10.5 deg. Yellow contours indicate the 3cr (0.15 Jy/beam) emission level 
of the ATLASGAL data. Red and blue ellipses show the Hii regions and SFCs, respectively. SLCs are shown with green filled diamonds. Similar 
maps for the Galactic plane between 10.5 deg < / < 21 deg and \b\ < 1 deg are shown in Appendix ICl 


- (19.2")^. This convolution was done for 

each region individually. At the distance of 5 kpc, the 19.2" reso¬ 
lution of the ATLASGAL translates to about 0.5 pc. We therefore 
do not resolve the dense cores ultima tely linked to star formation 
that have typically a size of ~0.1 pc (iMotte & HennebeH'3l2009l: 
iRagan et al.ll2013}) . 

When smoothing maps to a common distance, some of the 
smaller SLCs were washed out by strong emission gradients 
likely associated with nearby strong sources. This artificially in¬ 
creases the SLCs column densities. To minimize the effect, we 


inspected each SLC by eye, discarding those that were signifi¬ 
cantly affected by strong gradients. AppendixiDl shows the SLCs 
included in the final sample. 


The total number of regions studied in this paper, and the 
number of regions in each evolutionary class, are listed in Ta- 

blem 
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Table 1. Completeness of each evolutionary class 



Distance (kpc) 



Hii 

SFCs 

SLCs 

No class 

Total 

114 

184 

210 

107 

Known 

84 

126 

120 

— 

d < 5.5 kpc 

57 

99 

111 

— 

d > 5.5 kpc 

27 

27 

9 

— 

miss, d < 5.5 kpc^ 

20 

18 

102 

87 

miss, d > 5.5 kpc^ 

10 

8 

8 

20 

Studied 

57 

99 

3R 

— 


Notes. Only SLCs with KDA solved and, if more than one distance 
estimate, agreement between different literature sources. Number of 
regions lost due to lack of distance estimates. We assume homogeneous 
distribution of the sources along the Gala ctic plane area studied. We 
only studied isolated SLCs (see Sect. 12.21 ) 


Fig. 2. Distance distribution of the molecular cloud regions. Black solid 
line shows the total number of regions. Red dotted line shows the HII 
regions, the blue dashed line the SFCs, and green filled area the SLCs. 
Black dashed vertical line at 5 kpc shows the common distance to which 
we have smoothed the data. 



Fig. 3. Artist impression of face-on view of the Milky Way (R. Hurt, 
SSC-Caltech, MPIA graphic, Tackenberg et al. 2012). Hii regions are 
shown as red circles, star-forming clouds as blue circles and starless 
clumps as green circles. Circle sizes are proportional to region sizes. 
The right panel shows a zoom to the region enclosed by the black rect¬ 
angle in the left panel, where the source density is highest. 


2.3. Column density and mass estimation 

Column densities of molecular gas were calculated via 

AH2[cm“^] =--, 

B;t(TDust)fimuKQ' 


(3) 


where and B^(T) are respectively the flux and the blackbody 
radiation as a function of temperature, T, at 870yum. The quan¬ 
tity yu is the mean molecular weight (assumed to be 2.8) of the 
interstellar medium per hydrogen molecule, mu is the mass of 
the hydrogen atom, Q is the beain solid angle, and R = 154 is 
the gas-to-dust ratio (lDrainell201ll) . We used a dust absorption 


Table 2. Literature sources from which distances were obtained 


Reference 

Hii 

SFCs 

SLCs 

1 

7 

50 

11 

2 

19 

6 

18 

3 

23 

1 

7 

4 

39 

23 

18 

5 

7 

11 

5 

6 

14 

26 

5 

7 

— 

13 

— 

8 

17 

— 

5 


References, fl) lEllsworth-Bowers et al.l (120131 

): 12) iTackenberg et al.l 

J 2 OI 2 

): G) lUrauhart et al.l (120131): 14) 1 

IWienen et al.l ( 

20121): 

200^; 

(5) 

Roman-Duval et al.l (120091): (^6) 

Marshall et al.l ( 

('7)ISimon et al.l (120061): r8)IWalsh et al.l (11997 

). 


coefficient a: = 1.8 5 cm^g~^ at 870//m, wh i ch wa s calculated by 
interpolation of the lOssenkopf & Henningl (1 19941) dust model of 
grains with thin ice mantles and a mean den sity of ^ = 10^ cm~^ . 
We assumed T = 15 K for SLCs and SFCs (IWienen et al.ll2012h . 
in agreement with pr evious dust temperature estimations within 
infrared dark clouds (iPeretto & Fullejl201Ql IRDCs) and in en¬ 
velop es of star-forming cores (IStutz et al .1120101: iLaunhardt et al.l 
l2013l) . For Hii regions we assumed T = 25 K. This dust tem¬ 
perature is in agreement with the average dust temperatures 
in PDR regions surrounding H ii regions {T = 26 K), where 
most of the FIR-submm du st emission of these objects comes 
from (I Anderson et al.ll2QT^ . It also agrees with the mean tem¬ 
perature found in the central re gion of NGC6334 { T ~ 24 K), 
that is an expanding Hii region (iRusseil et al]l2013l) . For a bet¬ 
ter comparison with previous works, we present the column 
density data also in units of visual extinct ion using a conver - 
sion: A^H 2 = 0.94 x 10^^ Ay cm“^ mag“^ (iBohlin et al.lll978l) . 
The rms noise of the ATLASGAL data (50mJy) corresponds to 
Ay = 4.5 mag for both the SFCs and SLCs and 2.2 ma^ for 
Hii regions. No saturation problems were found in the ATLAS¬ 
GAL survey. The optical depth is « 1, therefore our measure¬ 
ments do not su ffer from optical dept h effects in the high-column 
density regime (ISchuller et al.ll200% . 


^ The difference in the rms values in terms of Ay is due to the temper¬ 
atures assumed for each evolutionary class. 
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We estimated the total gas mass of each region from dust 
continuum emission, assuming that emission is optically thin: 




Rd^F^ 

Bx{TDust)K 


(4) 


where d is the distance to the region. We assume the same val¬ 
ues for the other listed quantities as we assume for the col¬ 
umn density determination (Eq.O. Masses of the regions cover 
three orders of magnitude (Fig.|4l). The masses of the SLCs span 
0.2-4X10^ Mo, the SFCs 0.3 -15 x 10^ Mq, and the H ii regions 
0.2 - 200 X 10^ Mq. Farger masses for H ii regions and SFCs are 
expected since both have much larger sizes than SFCs (Fig.[T]). 

The derived mass and column density values depend on the 
assumed dust properties, specifically on Arsvoyum, R and Toust- 
Both ATsTOyum and R are subject to uncertai nties: K^ioum values 
differ by ~ 1 dex in different dust models (IShirlev et al.ll2005l 
I 2 OIII) . Eq.[3]and Eq.|4|assume isothermal clouds. This is clearly 
an oversimplification, increasing the uncertainty in the derived 
masses. Mass depends also on making uncertainties in dis¬ 
tance a major contributor to the absolute uncertainties. If we 
adopt Ad ~ 0.5 kpc, nearby regions will be more affected by 
distance uncertainties (50% at 1 kpc) than the most distant re¬ 
gions (10% at 5 kpc). This assumption agrees with th e distance 
uncertainties reported by (iRoman-Duval et alll2009l) . We note 
that the absolute uncertainty in our derived column densities is 
very large, potentially larger than a factor 10. The relative uncer¬ 
tainties between the evolutionary classes can be influenced by 
the different temperature assumptions or intrinsic differences in 
the dust properties. The isothermal assumption introduces dif¬ 
ferences in the low-column density regime of the WPDFs, but it 
has negligible effect in shaping the column density distribution 
at high-column densities (see App. IB. In the case of dust proper¬ 
ties, we have no knowledge about any observational-based study 
suggesting changes in them in molecular clouds at different evo¬ 
lutionary phases. We therefore assume that the dust properties do 
not introduce relative uncertainties between the three molecular 
cloud classes defined. 


2.4. Physical properties of the evolutionary classes 

We define and analyze in this work three distinctive evolution¬ 
ary classes of objects: SFCs, SFCs and Hit regions. The ob¬ 
jects in these classes are different in their physical characteris¬ 
tics. These differences originate dominantly from the fact that the 
HII regions and SFCs are typically extended regions (i.e., molec¬ 
ular clouds or even cloud complexes), while SFCs are smaller, 
"clump-like" structures. We quantify here the basic physical 
properties of the objects in our three evolutionary classes. The 
properties are also listed in Table (H 

Figure |4] shows the mass distribution of our regions. The 
mass distribution of H ii regions spans 3 dex from 10^-10^ Mq. 
SFCs have masses of 10^ - 10"^ Mq, and SFCs show the most 
narrow mass range, 10^ - IO^Mq. 

The spread of the distribution of mean column densities is 
Ay = 3 - 25 mag and it peaks at Ay = 7 mag (see Fig. [5]). The 
Ay distribution differs in each evolutionary class. While most 
SFCs and SFCs have Ay < 10 mag, a considerable number of 
Hii regions show Ay > 10 mag. We note that the mean column 
densities of our sample are overestimated due to the spatial fil¬ 
tering of ATFASGAF, and so is the peak of the Ay distribution. 
This effect is more important in H ii regions and SFCs since they 
have larger areas and hence larger fraction of diffuse material 
that is filtered out than the SFCs. 
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Fig. 4. Mass distribution of the molecular cloud regions. Filled green 
area shows mass for starless clumps, dashed blue line shows star¬ 
forming clouds and red dashed line shows H ii regions. Masses are given 
in solar mass units. 





Fig. 5. Mean column density. Ay, and size distribution of all the regions. 
In the scatter plot we show the H ii regions in red, the SFCs in blue and 
the SLCs in green. The histograms show the Ay and size distributions 
of each evolutionary class (same colors) and the whole sample (black). 


Figure [ 5 ] also shows the size distribution of each class and of 
the total sample. The SLCs have the smallest sizes of the sample 
with a mean size of 1.4 pc and a range of sizes between 1 pc and 
2.5 pc. The range of sizes of the SFCs is 2-15 pc, with a mean of 
5.3 pc. The Hii regions have the largest mean size of the three 
evolutionary classes, 7 pc, and also the largest spread, 2-18 pc. 
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Table 4. Mean physical properties of the evolutionary classes 



Hii 

SFCs 

SLCs 

Mass [Mef 

18x10^ 

2.7x10^ 

1.2x10^ 

Ay [mag] 

10.6 + 6.2 

7.0 ±2.5 

8.2 + 4.2 

Size [pc] 

7.0+ 4.0 

5.3 + 2.7 

1.4+ 0.6 


3. Results 

3.1. Column density distribution 

We use the column density data to study the column density dis¬ 
tributions of the regions. In the following, we first analyze the 
A^-PDFs and DGMFs. We then examine the relationship between 
the total mass and the column density distribution of the regions. 

3.1.1. /V-PDFs 

We first analyze the total WPDFs of the three evolutionary 
classes. To construct the WPDFs, used the mean normal¬ 
ized column densities v = In (Ay/Ay) (see Eq. [T]) of each re¬ 
gion. We calculated Ay as the mean column density of all the 
pixels of each region. The resulting WPDFs were then stacked 
together to form the total WPDFs, shown in Fig. [6] as the black 
histogram. The three classes show clearly different WPDFs: the 
HII regions have the widest (or shallowest) WPDF, followed by 
a slightly narrower (or steeper) WPDF of SFCs. The SLCs have 
the narrowest A^-PDF. 

Interpreting the low-column density shape of the A^-PDFs re¬ 
quires taking into account two issues. First, ideally the A^-PDF 
should not be affected by how exactly the field-of-view toward 
an individual region is cropped, i.e., it must include all column 
density values above a given level. Second, one must ascertain 
that the pixels are not dominated by noise or contamination from 
neighbouring regions. To fold these two limitations into one, we 
define a ’’reliability limit” of the A^-PDFs as the minimum col¬ 
umn density value above which all regions of the evolutionary 
class are well defined by a closed emission iso-contour (see Ap¬ 
pendix |D1|E| and 10. These levels are v = -1.5, -0.75, and 0 for 
HII regions, SFCs, and SLCs, respectively. These levels corre¬ 
spond typically to Ay = 2,4, and 9 mag, all at least l-cr above 
the noise level (50mJy; see Section (TS). The larger reliability 
limit in SLCs originates from the fact that they are embedded 
in HII regions and SFCs, being surrounded by emission levels 
higher than the map noise. The total number of pixels above 
these limits are 20 x 10^, 9 x 10"^, and 10"^ for H ii regions, SFCs, 
and SLCs, respectively. We note that this definition of the relia¬ 
bility limit is very conservative; it is set by the lowest iso-contour 
above which all regions of the evolutionary class show a closed 
contour. Most regions, however, have this limit at lower v val¬ 
ues. We also note that systematic uncertainties such as the dust 
opacity uncertainty do not affect (or, are unlikely to affect) the 
relative shapes of the three classes with respect to each others. 

To quantify the shapes of the A^-PDFs shown in Fig. [6j we 
fit them, equally sampled in the log space, with a combination 
of log-normal (see Eq. [T]) and power-law (p(s) oc cs^) func¬ 
tions. We used five free parameters in the fit: the width (cr^) and 
mean (ja) of the log-normal function, the slope (p) and constant 
(c) of the power-law, and the break point between both func¬ 
tions (st). Furthermore, molecular cloud masses should be re¬ 
covered when integrating the fitted function, representing an ex¬ 
tra boundary condition to the fit. The fitting range was defined 


as all V values larger than the reliability limit. We weighted the 
data points by their Poisson noise. We obtained the uncertainties 
of the fitted parameters by fitting the A^-PDFs using different bin 
sizes (ISadavov et al.ll2014l) . Results are summarized in Table [5l 

SLCs are well described by a log-normal A^-PDF (ct^ slc = 
0.5 + 0.1). Even though the peak of the A^-PDFs is below the 
reliability limit, it is well constrained by the fit because of the 
normalization factor in Eq. [T] The A^-PDFs of H ii regions and 
SFCs are inconsistent with a single log-normal function; they 
are better described by a combination of a log-normal function 
at low column densities and a power-law function at high col¬ 
umn densities. The low-column density log-normal portion of 
Hii regions is wider ((T 5 ,hii = 0.9 + 0.09) than that of SFCs 
(u~s,SFC = 0.5 + 0.05). The mean-normalized column densities 
at which the A^-PDFs transition from log-normal to power-law is 
similar in both classes, Hii regions and SFCs: St = 1.0 + 0.2. 
We also find differences in the power-law slopes of the N- 
PDFs. The power-law slope is clearly shallower for H ii regions 
(p = -2.1+ 0.1) than for SFCs (p = -3.8 + 0.3). 

3.1.2. Dense Gas Mass Fraction 

In Section [T] we defined the DCMFs as the fraction of gas mass 
enclosed by regions with M(Ay > Ay), relative to the total mass 
of the cloud (see Eq.O. Fig.|7]shows the mean DCMFs of each 
evolutionary class. Cenerally, H ii regions exhibit larger reser¬ 
voirs of high-column density gas than the SLC and SEC regions. 

We quantified the shapes of the mean DGMFs by fitting 
them with a combination of exponential (cx e^^^) and power- 
law (oc Ay) functions, leaving both exponents and the breaking 
point as free parameters and weighting each point by the Poisson 
standard deviation. Fit errors were calculated as in Section IT. 1.11 
resulting in parameter value uncertainties of 10%-15%. While 
the mean DGMF of the SLCs is well fitted by an exponential, 
HII regions and SFCs transition from an exponential to a power- 
law shape at Ay > 20 mag. This change is evidently linked to 
the change from log-normal to power-law shape in the A^-PDF 
because the DGMFs are an integral of the A^-PDF. H ii regions 
show the shallowest mean DGMF (a = -0.06), followed by 
SLCs (a = -0.11) and SFCs (a = -0.14). In the power-law por¬ 
tion of the DGMFs, Hii regions are also shallower (J3 = -1.0) 
than SFCs (J3 = -2.1). The amount of mass enclosed by the 
power-law DGMF is 30% of the total mass in H ii regions, al¬ 
most a factor of three lower, 10%, for the SFCs. 

The mean DGMF of Hii regions above Ay = 300 mag is 
dominated by regions #4, #55 and #122 (see Table IF.3I) . This 
fiat tail is built up by less than 1% of the pixels in each of the 
mentioned regions, and hence, is not representative of the whole 
Hii sample. 

3.1.3. Relationship between the region’s mass and column 
density distribution 

Does the dense gas mass fraction of a region depend on its mass 
an d from therein, affe ct the SFR - cloud mass relation presented 
by iLada et al.l (l2012h ? We analyze the DGMFs of each evolu¬ 
tionary class divided in five mass intervals (listed in Table [6]) 
that have at least 9 regions each to answer this question. 

Figure [8] shows the mean DGMFs of each mass interval for 
the three evolutionary classes. In all evolutionary classes, most 
massive regions have shallower DGMFs than those of less mas¬ 
sive regions. We fit the mean DGMFs with exponential and 
power-law functions as described in Sect. 13.1.21 Most DGMFs 
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s — ln(Av^/Av^) 



s — ln(Av^/Av) 




— <7 = 0.49 



S = \n{Ay/Ay) 


Fig. 6. Total mean normalized column density PDFs of Hii regions (left), SFCs (center) and SLCs (right). All panels: horizontal axis show mean 
normalized column densities, 5 = \n{Av/Av). Vertical error bars show Poisson standard deviation, cFpoisson ^ V^- The best-fit curves assuming 
a combination of log-normal and power-law functions are indicated, respectively, by red and green solid lines, with the fit errors indicated as 
shaded regions of the same colors. The gray shaded regions indicate data below the reliability limit. These data were excluded from the fit. Blue 
shaded regions show the range of values obtained for the mean-normalized column density value at which V-PDFs deviate from a log-normal to a 
power-law like function, Sf 





Fig. 7. Mean DGMFs of Hii regions (left), SFCs (center) and SLCs (right). Solid colored lines show mean normalized DGMFs. DGMFs were 
normalized to the reliability limit of each evolutionary class: Ay = 2, 4, 9 mag for Hii regions, SFCs and SLCs, respectively. Colored dashed 
lines show the fit of the DGMFs with exponential functions. Dashed-dotted colored lines show fits with power-law tails in the higher Ay range. 
Grey shaded regions show statistical poisson errors of the DGMFs. Small box in left panel shows the whole mean DGMF of Hii regions up to 
Ay = 1000 mag. 





Fig. 8. Mass-binned average DGMFs for each evolutionary class. Each line shows the DGMF for each of the mass intervals listed in the corre¬ 
sponding panel and defined in Table Dotted lines, dotted-dashed lines, dashed lines, and so lid lines progress from less to most massive bins, 
respectively. The DGMFs were normalized following the procedure described in Section [3. 1.2 1 and shown in Fig.[7l 


could not be fitted well with the combination of both functions 
over their entire column density range. Only the DGMFs of the 
most massive SFCs and H ii regions required two component 
functions; DGMFs of less massive regions are well described 
by an exponential alone. Exponents derived from this analysis 
are shown in Table [71 

In all evolutionary classes, the exponent of the exponential 
function, a, increases with mass (see TableO. In order to further 
investigate this correlation, we repeated the same fitting proce¬ 
dure for each individual region. Results are shown in Fig.O Top 


panel of Fig. shows the relationship between a and the mass 
of each region: a oc ^-0-43+o.o5^ a correlation coefficient 

r = 0.64 and a significance value p = 0.18. The fit parameters 
and their errors were obtained from a Monte-Carlo simulation 
of 10^ cycles. On each cycle we selected a random sample of 
points and fitted the resulting data set. We adopt the average of 
the best fit parameters obtained on each cycle as the best fit val¬ 
ues and their standard deviation as the error of the fit. It could 
be argued that the correlation between a and mass is dominated 
by the most massive (M > 3 x 10"^M©) Hii regions. To estab- 
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Table 5. Results of the best-fit parameters to the total V-PDFs and 
DGMFs. 



A^-PDFs 

DGMFs 


(T/ 


sp 



Hit 

0.9+0.09 

-2.1+0.1 

1.0+0.2 

-0.06^> 

-1.0 

SFCs 

0.5+0.05 

-3.8+0.3 

1.0+0.2 

-0.14 

-2.1 

SLCs 

0.5+0.1 

— 

— 

-0.11 

— 


Notes. Standard deviation of the log-normal portion of the V-PDFs. 

Slope of the power-law portion of the V-PDFs. Transition from 
log-normal to power-law portion of the V-PDFs in mean-normalized 
column densities. Slope of the exponential portion of the DGMFs. 

Slope of the power-law portion of the DGMFs. ^ Relative errors of 
DGMFs account for 10%. 

Table 6. Mass intervals of each evolutionary class in Fig. [8] 


[Mo] 

Hii 

SFCs 

SLC 

M,-,i 

— 


<lo^ 

Mi. 2 

1 - 2x 10^ 

1 - 2x 10^ 

1 - 2 X 10^ 

Mi. 3 

2 - 5 X 10^ 

2 - 5 X 10^ 

2 - 5 X 10^ 

Mi. 4 

0.5 - 1 X 10^ 

0.5 - 1 X 10^ 

— 

Mi. 5 

o 

A 

— 

— 


Table 7. Slopes of the exponential and power-law fits to DGMFs, a, p 
for the mass ranges presented in Tabled 




Mi. 2 

Mi. 3 

Mi. 4 

Mi. 4 

a(Hii ) 

— 

-0.25 

-0.22 

-0.06 

-0.04 

a(SFCs) 

-0.29 

-0.20 

-0.18 

-0.10 

— 

a(SLC) 

-0.32 

-0.19 

-0.09 

— 

— 

AH II ) 

— 

— 

— 

-0.59 

-0.61 

ASFCs) 

— 

— 

-1.03 

— 

— 


lish whether the correlation strength depends strongly on these 
few massive clouds we also explored this correlation without 
those extreme points. The correlation coefficient is somewhat 
lower in this case r = 0.56 and a significant value p = 0.21. 
However, there is no significant difference in the resulting fit 
(a oc ^-0-40+0.08^ Power-law exponents of DGMFs also exhibit 
a correlation with mass: p oc M 0-16+0.03 middle panel of 
Fig. O. The larger scatter seen in the data from the exponential 
fits relative to that seen in power-law fits may indicate that the 
power-law regimes of DGMFs are much better constrained than 
the exponential regimes. 


4. Discussion 

4.1. N-PDFs as a measure of the evolutionary stage of 
objects. 

The total A^-PDFs of different evolutionary classes exhibit clear 
differences; these differences can be linked to differences in the 
mechanisms that drive the evolution of objects within the var¬ 
ious classes. The A^-PDF of SLCs is well described by a sin¬ 
gle log-normal function (see Fig. This agrees with previ¬ 
ous o bservations of starless low-mass clouds (iKainulainen et al. 
2009 ) or star le ss regions of s t ar-form ing clouds (ISchneider et al. 
2012 . 120131: iRusseil et~^ l2013l) . In particular, this sim¬ 
ple log-normal form agrees with predictions for turbulence- 





Fig. 9. Top: relationship between the mass of the regions [M©] and the 
exponent of the exponential fit to the DGMFs (aexp)- Black solid line 
shows the best fit to the data and the shaded region shows its cr error. 
The dotted and dashed lines show the best fit when the most massive 
HII regions are removed. Colors indicate the evolutionary class of each 
point as indicated. Middle: relationship between mass of the regions 
and the slope of the power-law range of the DGMFs. Black line shows 
the best fit to the data and the shaded region shows its cr error. Bottom: 
relationship between the mean gas mass surface density of the MCs, 
Emass. and the dense gas mass fraction of gas, /dg = The 

crosses show the /dg obtained integrating the exponential regime of the 
DGMF s of e ach region in the range Ay = 0-7 mag (see third paragraph 
in Sect. 14.21 ). Black line shows a linear fit to the data in the ra nge Emass = 
50-200 M© pc~^. Vertical dashed line at Emass = 116 Mopc"^ jLada et al.l 
1201 011201% indicates the threshold for the dense gas. 


dominated media from numerical simu lations (iPadoan et al.l 
Il997l:lvizquez-Semadeni & Garcfall200ll) . 
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In contrast, the total A^-PDFs of star-forming clouds, i.e., 
H II regions and SFCs, show two components that can be de¬ 
scribed by log-normal and power-law functions. The power-law 
components of the A^-PDFs of the H ii regions are shallower than 
those of the SFCs. Previous studies have found that the power- 
law slopes are within p = [-1.5,-3.3], with shallower slopes re¬ 
lated to most active star-fori ning regions. In the o nly study with 
a resolution similar to ours, iRusseil et al.l (12013h found a non- 
star-forming region in NGC 6334 to have a steep A^-PDF slop43 
(p = -5.7), moderately star-forming regions to have shallower 
slopes (p = -3.3, -3.0), and an Hii region to have the shallow¬ 
est slope (p = -1.5). Their trend to have shallower A^-PDF in 
a cloud region that contains an H ii region is similar to what we 
find in our work. 


Theories and simulations that consider turbulent gas under 
the influence of gravity predict power-law-like tail s for A^-PDFs 
with e xponents comparable to wh at is observed (iKritsuk et al.l 
1201 iL iFederrath & Klesseiil l2013l) . possibly featuring batten¬ 
ing of the power-law over time-scales relevant for star forma¬ 


tion (iBallesteros-Paredes et all201ll:lFederrath & Klessenll2013l: 


iGirichidis et al.ll2014l) . iKritsuk et alJ (1201 ll) showed that a col¬ 
lapsing spherical cloud with a power-law density distribution, 
p oc r“^, will have a power-law A/^-PDF with a slope of p = 
-2/{k - 1). The power-law slopes that we observe (see Table [5]) 
indicate k = 1.9 and 1.5 for the Hii regions and SFCs, respec¬ 
tively. The former is very close to the va lue a: = 2 of a col¬ 
lapsing isothermal sphere (IShu et al.ll 19871) . suggesting that the 
density distribution of H ii regions may be dominated by self¬ 
gravity. The value of a: = 1.5 we find for SFCs can also be 
indic ative of a collapse slow ed down by turbulence-driving ef¬ 
fects (iGirichidis et al.|[2QT^ . We note a caveat in this analysis. 
Our SFCs and Hii regions are unlikely to be close to spheres 
and their large sizes make them unlikely to be under general 
free-fall collapse. However, these regions are composed of nu¬ 
merous smaller ATLASGAL clumps (see Fig. [T]) that may be 
closer to spherical symmetry, and we are averaging the emission 
of all these smaller clumps. Indeed, the density profile expon ent 
of our SFCs is similar to that found by iBeuther et al.l (l2002h in 
a sample of small high mass star-forming objects, which corre¬ 
spond to our definition of SFCs. 


Recent works based on Herschel observations have 
explored possible effects of other processes (e.g. ionis¬ 
ing radiatio n or shock compression) on the A-PDFs o f 
Hn regions (ISchneider et al.ll2012L 120131: iTremblin et al.ll2014l) . 
iTremblin et al.l (12014l) report A-PDFs with two log-normal com¬ 
ponents. They relate the log-normal component at low column- 
densities to the turbulent motions of the gas and the component 
at high-column densities to ionization pressure. They also sug¬ 
gest that the presence of these double-peaked A/^-PDFs depends 
on the relative importance of ionizing and turbulent pressures. 
The total A-PDFs of our H ii regions, composed of 60 individ¬ 
ual regions, does not exhibit such behavior. This could originate 
from a combination of seve ral factors: i) the low -column den¬ 
sity component detected by ITremblin et al.l (l2014l) is at column 
densities of Ay < 6 mag. These column densities are generally 
filtered out from the ATLA SGAL data; //) t he siz e-scales of the 
molecular clouds studied in ITremblin et al.l (l2014l) and this work 
are different and it may happen that the ionisation front of the 
H II regions is not spatially resolved in our observations. 


^ iRusseil et all (1201 3h quote the equivalent radial density profile (a:), 
which can be related with the slope of the power-law tail of the A-PDF 
via p - -2/(a: - 1). 


The above models offer an attractive possibility to link 
the observed A-PDF tails to self-gravitating gas in molecu¬ 
lar clouds. However, it has not yet been shown observation- 
ally that the power-law parts would be definitely caused by 
sel f-gravity; an alter n ative i nterpretation has been proposed 
by iKainulainen et al.l (1201 lah who suggested that the overall 
pressure conditions in the clouds may play a role in produc¬ 
ing the observed power-law-like behavior in low-mass molecular 
clouds. 


4.2. Dense Gas Mass Fraction in molecular clouds 

With our cloud sample, we are able to study the DGMFs of 
molecular clouds over a relatively wide dynamic range of col¬ 
umn densities and separately in various evolutionary classes. The 
continuous DGMF functions (Eq. [2]) allow for a more complete 
census of the dense gas in the clouds than the analysis of the ra¬ 
tios of two tracers, e.g., of CO emission and dust emission. We 
find that the DGMFs of H ii regions are shallower than those of 
SFCs and SLCs. This suggests a direct relation between the star¬ 
forming activity of molecular clouds and their relative dense gas 
mass fraction. S imilar results have also been previously found in 
nearby regions dEada et al.ll200^l20 id: IKainulainen et al.ll2009l: 

Kainu^lainen & TanI 1201 3l) and filmentary clouds dAndre et al.l 

2f)KilK 

We detect a clear correlation betw een the DGMF slope and 
cloud mass (see Fig. [8]). Previously, iBattisti & Heverl (1201 4l) 
found no correlation between molecular cloud mass and the 
dense gas fraction in a large sample of molecular clouds. They 
defined the dense gas fraction as the ratio of dust emission- 
derived mass, traced with 1 mm flux, above Ay = 9.5 mag 
to CO-derived mass: /dg = Their result imply 

that there is no correlation between the mass of CO-traced gas 
(Ay ~ 3 - 8 mag) and the mass of gas at column densities 
Ay > 9 mag. Unfortunately, we do not measure the CO mass of 
our MCs and therefore we cannot dir ectly compare our results 
with those of IBattisti & HeveJ (l2014t) . The correlation we find 
between the molecular cloud masses and the slope of DGMFs 
suggests that the dense gas fraction depends on the mass of mod¬ 
erately dense gas (Ay > 10 mag) rather than the CO mass of the 
clouds. 

iLada et aP (120101 l2012h suggested that star formation rates 
depend linearly on the amount of dense gas in molecular clouds: 
^SFR oc /oG^mass, wit h fpG = M(Av > 7 ma g)/Mtot. Combin¬ 
ing this relation with iGutermuth et al.l (1201 ll) . who derived the 
relation Esfr oc suggests /dg oc Emass- We find that this 

correlation indeed exists in the range Smass = 50 - 200Mopc“^ 
(see Fig.O. At higher surface densities the relationship flattens 
at /dg = 0 .8, suggesting that the maximum amount of dense 
gas that a MC can harbor is around 80% of its total mass. Con¬ 
sequently, the maximum Esfr of a molecular cloud is reached 
at /dg =0.8. This value depends on the definition of the col¬ 
umn density threshold (A*) of the dense gas becoming lower for 
higher values of A*. The spatial filtering of ATLASGAL data 
(see Appendix|A]) results in overestimated /dg values. We there¬ 
fore propose /dg = 0.8 as an upper limit to the actual maximum 
/dg of a MC. The overestimation of the /dg values derived above 
can be studied using DGMFs. In general, DGMFs have been 
shown to follow a n exponential function, o c down to low 
colum n densities (IKainulainen et al.l 120091: IKainulainen & TanI 
l2013l) . We adopted the a values calculated in Section [7.1.31 and 
integrated the exponential DGMF in the range Ay = 0-7 mag 
to obtain an estimate of /dg- The result is shown with crosses 
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Table 8. Statistics of /dg in this work and in jLada et al.ll2012l) . 


Aftot 

/dg o-K/dg) # of regions 

This work 



< 0.8 X 10^ 

0.24 

0.22 

13 

0.8 - 2.2 X 10^ 

0.29 

0.73 

42 

2.2-6.0 X 10^ 

0.31 

0.64 

39 

6.0 - 17 X 10^ 

0.41 

0.50 

27 

17-46X 10^ 

0.58 

0.33 

10 

> 46 X 10^ 

0.83 

0.16 

4 

Entire range 

0.39 

0.71 

135 

(Lada et al. 2012) 

0.8- 100 X 10" 

0.11 

0.56 

11 


in bottom panel of Fig [9j The mean overestimation of /dg in 
SFCs and Hit regions is 2 and -1.3 respectively. We did not in¬ 
clude the SLCs in this experiment because their reliability limit 
is Ay = 9 mag and therefore they have /dg = 1 (i-e. all its mass 
is enclosed in regions Ay > 7 mag). 

Our data can also help to understand the S FR - dense gas 
mass relation suggested bv iLada et al.l (l2012h . The S FR - dense 
gas mass relation shows significant scatter of star forn iation rates 
for a g iven dense gas mass, about 0.6 dex (see Fig. 2 lLada et^ 
(l2012h and Fig.[T0|). This scatter shows that not all clouds with 
the same amount of dense gas form stars with the same rate. To 
gain insight into this, we calculated the dense gas ma ss fractions 
and sta r formation rates for our regions as defined bv iLada et al.l 
(I2012h . i.e., SFR = 4.6 x lO-^focMtot M suny^ ^ • Figure [TOl 
show s the SFR - den se gas mass relation with data points 
from iLada et al.l (l201Ql) . The figure also shows the mean SFR 
of our regions in six mass bins, with error bars showing the rel¬ 
ative standard deviation of /dg- The standard deviations are also 
listed in Table [8j The relative standard deviation of /dg over the 
entire mass range of our regions is 0.71, which is s lightly higher 
than the relative scatter of SFR in lLada et~^ (l2012h . /dg = 0-56. 
We conclude that the scatter in star formation rates for a given 
dense gas mass can originate from differences in dense gas frac¬ 
tions, i.e., in the total masses of clouds for a given dense gas 
mass. This, in turn, suggests that the dense gas mass is not the 
only ingredient affecting the star formation rate, but the lower- 
density envelope of the cloud also plays a significant role. How¬ 
ever, we note the caveat that ATLAS GAL filters out low-column 
densities , which may ma ke the dense gas fra ctions we derive not 
comparable with those in iLada et al.l (l2012h . derived using dust 
extinction data. 

4.3. Evolutionary time-scales of the evolutionary classes as 
indicated by their N-PDFs 

If A^-PDFs evolve during the lives of molecular clouds, could 
they give us information about the evolutiona ry timescales of 
the clo uds in the three classes we have defined? iGirichidis et al.l 
(l2014l) have developed an analytical model which predicts the 
evolution of the p-PDFs of a system in free-fall collapse. They 
estimate the relative evolution time-scale, from the free-fall 
time at the mean density, p, of the molecular cloud, %(p), and the 
density at which the p-PDFs begin to show a power-law shape, 
Ptail 


tE = (5) 

V Ptail 



Solid black vertical lines show the standard deviation, cr, for each mass 
bin for our study. Black dashed line shows the constant value /dg = L 


The model also predicts the mass fraction of gas in regions with 
P > Ptail- We denote this mass as Mdense- 

Since our work is based on column densities instead of vol¬ 
ume densities, we need to write Eq. [5]in terms of column den¬ 
sities. To this aim, we assume a ratio between 2D and 3D vari¬ 
ances, R = (^^ni<n>/^'Iip' relation is also valid for the ra¬ 
tios p/ptaii and Ay/Ay^^ = where A^^^^ is the column density 
value at which the WPDF becomes a power-law and St is the 
mean normalized A^^^^ Then, Eq.[5]can be written as 




f O.2 

Vr 


e ^nfi{Ay). 


( 6 ) 


This equation allows us to estimate the evolutionary time- 
scale of a molecular cloud using two observable quantities, 
namely Ay and A^ y \ The factor R is still not well constrained. 
iKainulainen et al.l (12014l) obtaine d observationally a value of 
R - 0.4 while iBrunt et al.l (1201 oh obtained/? = [0.03,0.15] in 
MHD turbulence simulations without gravity. In the following 
we use the observationally derived value, R = 0.4, to estimate 
the time-scales of our three evolutionary classes. We estimate 
the uncertainty in the time-scales as the relative uncertainty be¬ 
tween R = 0.4 and /? = 0.15, which is roughly 30%. 

We find that our H ii and SFCs classes have evolutionary 
time-scales = ^e,sfc = 0-4 ± O.Hff and their relative mass 
of gas in regions with s > St are Mdense,Hii ~ 30 + 0.05% and 
Afdense,SEC ~ 10 ± 0.06% where the uncertainties were obtained 
from the l-cr uncertainties in St (see Section [3.1.11) . Since the N- 
PDF of the SLCs has no power-law tail, we calculated an upper 
limit of their evolutionary time-scale by using the largest extinc¬ 
tion in their A^-PDFs as a lower limit, St > 1.2. We obtained 
tE,SLC < 0.3 ± O.Hff. 

The above relative time-scales can be used to estimate abso¬ 
lute time-scales if the free-fall time is known. We es timate the 
free-fall time of each evolutionary class as = ^SnlSlGp. 
The mean density of each class was estimated using their 
mean masses and effective radifl and assuming spherical sym- 


^ We define the effective radius as the radius of a circle with the same 
area as the projected area of a given cloud. 
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Table 9. Evolutionary time-scales 



tE [tff] 

tE [Myr] 

4^dense [^] 

p [cm 

Hii 

0.4+ 0.1 

0.7 + 0.2 


0.3 X 10^ 

SECs 

0.4 ±0.1 

0.3+ 0.1 

1 G+0106 

^ -0.04 

1.5 X 10^ 

SLCs 

0.3+ 0.1 

<0.1 + 0.03 

— 

4.7 X 10^ 


metry. We find that the mean evolutionary time-scale for our 
H It regions is tE,uii = 0.1 ± 0.2 Myr, and the time-scale of SFCs 
is tE,SFC = 0-3 ± 0.1 Myr. SLCs have the shortest time-scales, 
tE,SLC < 0.1 + 0.03 Myr. We note that the absolute time-scales 
are measured using the onset of the gravitational collapse in the 
molecular cloud as r = 0 and that they were specifically esti¬ 
mated independently for each of the three classes of clouds de¬ 
fined in this work. 

Do the above results agree with previous time-scale esti¬ 
mations? The evolutionary time-scale of our SLC sample is 
within the range of collapse life-ti n ies de rived in other stud¬ 
ies. For example, 
of 6 X 10^ yr and 


Tackenberg et al.l (2Q12h derived a life-time 
Ragan et^ (12013), 7 - 17 x 10"^ yr for the 


starless core phase. Furthermore. ICsengeri et al.l (12014l) found 
a time-scale of 7.5 + 2.5 x 10"^ yr for massive starless clumps 
in the Galaxy using ATLAS GAL data. In all these studies, as 
well as in the present paper, the starless clumps are massive 
enough to be able to harbor high-mass star-forming activity. 
Similar evolutionary time-scales have also been found in re¬ 
gions th at are more likely to only form low-mass stars, e.g., in 
Perseus (IWalker-LaFollette et al.fl2014l) . The SFC evolutionary 
time-scale is close to recent age estimates of Class 0-rl proto¬ 
stars, ~ 0.5-0.4Mvr (iDunhani et al.ll20l4 Table 1). Hii regions 
are subject to other physical processes apart from gravity, such 
as Rayleigh-Taylor (RT) instabilities involved in the expansion 
of the HII regions and shocks due to stellar feedback. These pro¬ 
cesses make this simple evolutionary model hardly applicable to 
them and we therefore do not discuss the evolutionary time ob¬ 
tained for HII regions further. 

Finally, we mention several caveats associated with the time- 
scales derived above. The mean column density used in Eq. [6] 
corresponds only to the mean observed column density and not 
necessarily to the actual mean column density that should be 
used in Eq. [b] In addition, the factor, R, relating 2D and 3D 
variances of mean normalized densities is still not well con¬ 
strained. Eurthermore, this model assumes a single cloud under¬ 
going free-fall collapse. While this assumption can be true for 
the SLCs, it is unlikely to be the case in SECs. As mentioned 
in Section 031 we assume that the smaller ATLASGAL clumps 
which compose each SEC region are close to spherical symme¬ 
try. With these caveats, we only aim to study the evolutionary 
time-scales in terms of orders of magnitude. Considering these 
caveats, we conclude that the method of estimating evolutionary 
time-scale presented agrees with independently derived typical 
ages for SLCs and SECs. 


5. Conclusions 

We have used ATLASGAL 870 pm dust continuum data to study 
the column density distribution of 330 molecular clouds molec¬ 
ular clouds that we divide in three evolutionary classes: starless 
clumps (SLCs), star-forming clouds (SECs), and Hii regions. 
Our large sample of molecular clouds allows us to study their 
column density distributions at Galactic scale for the first time. 
We study the column density distributions of the clouds over a 


wide dynamic range Ay ~ 3 - 1000 mag, spanning a wide range 
of cloud masses (10^-10^ M©). In the following we summarize 
the main results obtained. 

- The total A-PDEs of SLCs is well described by a log-normal 
function with a width of about ~ 0.5. The total A-PDE of 
SECs and H ii regions show power-law tails at high column 
densities, with H ii regions having a shallower slope. These 
observations agree with a picture in which the density distri¬ 
bution of SLCs is dominated by turbulent motions. The SECs 
are significantly affected by gravity, although turbulence may 
still play a role in structuring the clouds. The density distri¬ 
butions of HII regions are consistent with gravity-dominated 
media. Our statistical sample shows that this picture, ear¬ 
lier observed in clouds of the Solar neighborhood, is relevant 
also at Galactic scale. 

- DGMEs of SLCs are well described by exponential functions 
with exponent aexp = -0.1. The DGMEs of H ii regions and 
SECs are better described by power-laws with exponents of 
13 = -1.0 and (3 = -2.1 respectively. The DGME shape de¬ 
pends on cloud mass, being shallower for the most massive 
clouds and steeper for the less massive clouds. This depen¬ 
dence exists in all evolutionary classes. 

- We find an approximately linear correlation /dg ^ ^mass for 
^mass = 50 - 2OOM0pc“^, valid for all evolutionary classes. 
This relation fiatens at /dg = 0.8 in MCs, suggesting that the 
maximum star-forming activity in MCs is reached at /dg = 
0.8. We also find that the intrinsic scatter of /dg is (-0.7 dex) 
is simila r to the scatter seen in the relation SER - dense gas 
mass of (iLada et al.ll2010ll2012h . This suggests that both, the 
dense gas mass and the lower-density envelope of the cloud, 
play a significant role in affecting the star formation rate. 

- We estimate the evolutionary time-scales of our three classes 
using an analytical model which pred icts the evolution of th e 
PDE of a cloud in free-fall collapse (ICirichidis et al.ll20T^ . 
We found < 0.1 Myr, tE - 0.3 Myr, and tE < 0.7 Myr 
for SLCs, SECs, and H ii regions, respectively. Both time- 
scales agree with previous, independent age estimates of cor¬ 
responding objects, suggesting that molecular cloud evolu¬ 
tion may indeed be imprinted into the observable A-PDE 
functions. Hii regions show a complexity of physical pro¬ 
cesses that make this model hard to apply to them. 
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Appendix A: Comparing ATLASGAL and Herschei 


Every observational technique to estimate A^-PDFs has its 
own limitations. The ATLASGAL data reduction process fil¬ 
ters out e xtended emiss i on fr om the maps in scales larger 
than 2.5' dSchuller et al.l 120091) The FIR emission observed 
by Herschei ( Pilbratt et al.l 1201 Cil) is very likely to be contam- 
inat ed by emission from dust unrelated to the cloud of inter¬ 
est (ISchneider et al.ll2015l) . We explore now how the WPDFs de¬ 
rived with Herschei and ATLASGAL differ. We do this for one 
example object of each evolutionary class using the Hii region 
M17 (#248), the SFC IRDC G11.11-0.12 (#54) and the SLC 
(#54c). 


The Herschei data of the SFC the SLC were taken as 
part of the Herschei guara nteed time key program Earliest 
Phase s of Star formation (iHenning et al.l l20ld: iRagan et ~sl\ 
12012L EPOS). The data of M17 was obtained from the Her¬ 
schei progra m Hi-GAL (iMolinari et al.l |2Q1Q|) . We used the 
three SPIRE (iGriffin et~^|2QlQ|) wavelengths (250 jtm 350 j/m 
and 500 jum), reduced u sing scanamorphos v23 (iRoussell 
l2013l) . and PACSJ_60yum (iPoglitsch et al.ll201^ , reduced using 
HIPE vl2 (lQttll2010l) . In the flux calibration process, the Planck 
zero-point correction was applied only to the Ml7 data. 

We derived the column density and temperature maps for 
each of the three selected regions through a pixel-to-pixel mod¬ 
ified greybody fit to the four Herschei continuum maps, all of 
them smoothed to a resolution of 36". For consistency with the 
ATLASGAL data analysis, we adopted the dus t opacity by in¬ 
terpolation of the lOssenkopf & Henningl (1 19941) dust model of 
grains with thin ice mantles and a mean density of ^ = 10^ cm“^. 
The mean uncertainty obtained in our greybody fitting technique 
for the temperature maps is ~ 2.5 K. The relative uncertainty of 
the column density maps is ~ 40%. The column density maps 
are shown in Fig. IA.il 

The area over which a molecular cloud shows significant 
emission is different in ATLASGAL and Herschei column den¬ 
sity maps. We face this issue by comparing the WPDFs over two 
difiTerent areas: the area over which ATLASGAL shows signif¬ 
icant emission, which we will refer to as dense gas area (white 
contours in Fig. IA.il and panels in the mid row in the same fig¬ 
ure). We also compare the A^-PDFs derived from the entire areas 
shown in Fig. IA.il 

We now describe how the A/^-PDFs derived from Herschei 
and ATLASGAL data look like. In the SFC and the H ii region 
Herschel-&QviNQ& A^-PDFs show a clear log-normal and power- 
law combination. This combination is seen in both cases of area 
selection: dense gas area and whole map. The ATLASGAL- 
derived A^-PDFs of the SFC and the H ii region also have power- 
law tails at high column densities but they do not show a log¬ 
normal distribution at low column densities. In both cases, Her- 
5'c/z^/-derived A^-PDFs do not probe regions with Ay < 10 mag. 
Both, the ATLASGAL-derived and Herschel-dQxiNtdN of 
the SLC are unfortunately dominated by noise, making a com¬ 
parison impossible. At 36" of resolution the SLCs do not have 
enough pixels for an analysis. 

The absence of column densities Ay < 10 mag in the Her¬ 
schei data-set is very likely related to the line-of-sight contam¬ 
ination. To compare A^-PDFs of the datasets without this con¬ 
tamination, we subtracted the background emission from the 
Herschei column density maps. We estimated the magnitude of 
the line-of-sight contamination averaging the Herschel-dtmtd 
column densities inside the white boxes shown in the top row 
of Fig. lA.ll We found = 11.2 + 1.2 mag in the SFC and 
SLC regions and = 7.0 + 1.6mag in the Hit region. The 
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background-subtracted A^-PDFs are shown in the fourth and fifth 
rows of Fig. lA.ll for the whole map and the dense gas area. 
The background subtraction significantly widens the Herschel- 
derived A-PDFs in the low column density regime. It has, how¬ 
ever, very small effect in the high column-density regime, which 
is slightly flattened. 

To estimate the difference between the fitted parameters in 
the datasets we fitted the power-law tails of the A-PDFs, fol¬ 
lowing the same procedure as in Section 13.1.11 The fits were 
performed in the column density regimes Ay > 30 mag and 
Ay > 40 mag in the SFC and H ii region respectively. In all cases, 
the power-law portion of ATLAS GAL is somewhat shallower 
(Phii,ag = -1.2, psFC,AG = -2.0) than that obtained for Her- 
schel (pmi,u = -L6, psFC,u = -2.3). The power-law slopes ob- 
tai ned in this work for the SFC are shallower than those reported 
bv ISchneider et all (l2014l) . The difference is probably caused by 
the different column density ranges used to fit the power-law in 
the works. When the background component of Herschel is re¬ 
moved, the power-law tails flatten and become much more sim¬ 
ilar to those observed by ATLASGAL = -1.2, ~ 

-1.9). Using only the ATLASGAL emission area (middle row of 
Fig. lA.lb or the whole map (bottom row of Fig. lA.lb makes no 
significant difference in the slope of the power-law tails obtained. 
We conclude that the high-column density power-law parts of 
the A-PDFs are in good agreement between ATLASGAL and 
Herschel. The agreement is even better when the background 
contamination component of Herschel is removed. Note that a 
background correction to the Herschel column densities is usu¬ 
ally necessary, as the diffuse Galactic dust component is signif¬ 
icant at the Galactic plane. Therefore, one should consider the 
background subtracted A-PDF as a better estimate of the A-PDF 
of the cloud. 

We identify the absence of log-normal components in the 
ATLASGAL-derived A-PDFs as an effect associated to the spa¬ 
tial filtering in the data reduction process. This effect is sig¬ 
nificant in both H ii regions and SFCs, being less important in 
denser regions of molecular clouds where SLC objects lie. Spa¬ 
tial filtering is clearly seen at column densities Ay ~ 10-20 mag 
in Fig. lA.ll where the Herschel-dtmtd A-PDFs shows a clear 
excess compared to the ATLAS GAL-derived A-PDFs. Despite 
the significant differences shown by the A-PDFs derived at 
low column density regimes, the power-law tails at high col¬ 
umn densities are in good agreement, showing the ATLAS GAL- 
derived A-PDFs marginally flatter distributions than the Her- 
schel-dtmtd A-PDFs. Similar results are obtained when the 
ATLAS GAL-derived and Herschel-dmy^d DGMFs are com¬ 
pared. 


To study the temperature effects on the resulting A-PDFs 
we used the Herschel derived temperature distributions in the 
previous section to reconstruct the ATLASGAL column den¬ 
sity maps of the same t hree r egions. The results of this exper¬ 
iment are shown in Fig. IbH In the H ii region, the isothermal 
assumption underestimates the low column density regimes of 
the A-PDF, which remain practically unaffected at Ay > 40 mag. 
The isothermal A-PDF of the SFC overestimates the low col¬ 
umn density regime and remain similar to the A-PDF of the Her- 
5'c/z^/-derived temperature distribution at Ay = 10 - 90 mag. The 
isothermal A-PDF in the SLC is shifted to lower column densi¬ 
ties. 

The isothermal assumption is therefore valid in the high col¬ 
umn density regime (i.e. in the power-law tail) of the H ii region 
and the SFC examples shown here. We note that the relative 
temperature uncertainties are larger in the coldest regions (T~ 
12 - 15 K) of molecular clouds (i.e. in the densest regions) and 
these uncertainties could also result in the underestimate of the 
A-PDF observed at Ay > 90 mag in the SFC. Unfortunately, 
we cannot quantify the possible differences in the shape of the 
isothermal and the Herschel-doxiytd temperature distribution A- 
PDFs of the SLC. The isothermal A-PDF therefore offers a more 
accurate reproduction of the Herschel-dtmtd temperature dis¬ 
tribution A-PDF in the column density regime of the power-law 
tail than at low-column densities. 

Appendix C: MIPSGAL 24)um maps with ATLASGAL 
contours and regions 

Please find this appendix in the journal version 

Appendix D: ATLASGAL maps from stariess 
ciumps 

Please find this appendix in the journal version 

Appendix E: Hii regions 

Please find this appendix in the journal version 

Appendix F: star-forming ciouds 

Please find this appendix in the journal version 


Appendix B: Effects of the isothermai assumption 
on the N-PDFs 

Obtaining column densities via dust emission maps at sub-mm 
wavelengths requires the use of the dust temperature (see Eq. a. 
When only one wavelength is available, as in the case of this 
paper, the most simple assumption is that the dust is isothermal. 
However, molecular clouds are not isothermal and the isothermal 
assumption can therefore generate artificial features in the col¬ 
umn density distributions of the maps derived with this method. 
When several wavelengths are available, as in the case of Her¬ 
schel observations, the line-of-sight averaged temperature and 
column density distributions can be simultaneously obtained via 
modified blackbody fitting to the FIR/sub-mm spectral energy 
distribution. 
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Fig. A.l. Top row: Herschel-dmvQd column density maps of M17 (Hii region), Gil (SFC) and #53c (SLC), in units of Ay. The white contours 
show the dense gas area (Ay >2 mag, 4.5 mag and 9 mag for the H ii region, SFC and SLC respectively). The white dashed boxes show the regions 
where the background contamination of Herschel has been calculated. Second row: V-PDFs as seen by Herschel (blue) and ATLASGAL (black) 
in the maps shown in the top row. The vertical error bars show the Poison standard deviation. The solid lines show the best fit to the power-law tail. 
Third row: ATLASGAL-derived and Herschel-d^nvo^d A-PDFs in the dense gas area. Fourth row: Background corrected ATLASGAL-derived 
and Herschel-dmy^d A-PDFs in the whole map area of the top row. The background emission was estimated as the mean column density in the 
dashed boxes of the first row, seen by Herschel. Bottom row: Background corrected A-PDFs evaluated in the dense gas area. 
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Fig. B.l. Isothermal ATLASGAL-derived A/^-PDFs (black) and A-PDFs derived using ATLASGAL emission maps together with Herschel-derived 
temperature maps (red). From left to right, Hii region, SFC and SLC. The vertical error bars show the Poison standard deviation. The solid lines 
show the power-law fit to the data in the column density range covered by the lines. 
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Table R3. Regions studied in this paper 


ID 

Type 

RA (J2000) 

DE (J2000) 

Radius (") 

Maj. Ax (") 

Min. Ax (") 

Angle (°) 

Distance (kpc) 

Dist. dispersion (kpc) 

References 

2 

SFC 

18:05:35.91 

-20:52:28.00 

90 

— 

— 

— 

3.2 

— 

1 

3 

SFC 

18:05:47.64 

-21:00:32.00 

90 

— 

— 

— 

4 

— 

1 

4 

HII 

18:06:15.12 

-20:31:36.50 

150 

— 

— 

— 

5.5 

0.4 

3 

5 

SFC 

18:06:46.98 

-20:59:04.10 

180 

— 

— 

— 

4.5 

— 

1 

7 

HII 

18:06:52.03 

-21:04:27.35 

180 

— 

— 

— 

4.8 

<0.1 

4,6 

8 

SFC 

18:06:53.75 

-21:18:46.50 

150 

— 

— 

— 

4.4 

— 

1 

11 

SFC 

18:07:34.17 

-20:26:03.30 

150 

— 

— 

— 

2.6 

— 

1 

13 

SFC 

18:07:42.18 

-21:23:02.20 

200 

— 

— 

— 

4.7 

0.3 

1 

16 

SFC 

18:07:55.85 

-20:28:26.60 

220 

— 

— 

— 

2.5 

— 

1 

27 

HII 

18:08:53.62 

-18:16:08.10 

350 

— 

— 

— 

3 

0.1 

4,2 

28 

SFC 

18:08:59.71 

-20:11:29.40 

150 

— 

— 

— 

2 

— 

4 

29 

HII 

18:09:23.23 

-20:08:38.60 

150 

— 

— 

— 

3.6 

0.1 

4,3 

31 

SFC 

18:09:24.38 

-20:01:35.70 

220 

— 

— 

— 

2.3 

0.4 

1,4 

30 

HII 

18:09:01.33 

-19:48:33.30 

— 

290 

250 

-60 

5.2 

— 

2 

33 

HII 

18:09:05.33 

-19:28:00.80 

150 

— 

— 

— 

2.5 

— 

3 

36 

HII 

18:09:37.96 

-20:19:16.10 

540 

— 

— 

— 

1.9 

0.2 

4 

37 

SFC 

18:09:06.10 

-21:03:37.20 

420 

— 

— 

— 

3.7 

0.3 

4 

38 

SFC 

18:09:22.35 

-21:15:35.70 

300 

— 

— 

— 

5.2 

— 

6 

44 

HII 

18:09:29.30 

-19:16:19.30 

— 

620 

500 

-60 

3.4 

— 

8 

47 

SFC 

18:09:54.44 

_19:44:46.90 

— 

450 

360 

-10 

3.5 

0.2 

1,4,2 

49 

SFC 

18:10:02.27 

-18:50:12.00 

280 

— 

— 

— 

3.4 

— 

1 

50 

HII 

18:10:05.43 

-20:59:11.90 

360 

— 

— 

— 

3.7^ 

— 

4,8 

51 

HII 

18:10:12.40 

-20:46:22.90 

— 

800 

260 

-50 

3.7 

0.2 

4,8 

54 

SFC 

18:10:26.70 

-19:20:56.20 

— 

920 

380 

-65 

3.5 

0.2 

1,4,3,2,6 

55 

HII 

18:10:54.13 

-19:52:35.30 

— 

860 

630 

-30 

5.2 

— 

6 

60 

HII 

18:10:51.72 

-17:55:57.40 

240 

— 

— 

— 

2.4 

0.3 

1,4,6 

61 

SFC 

18:10:54.72 

-20:32:50.30 

180 

— 

— 

— 

5.4 

— 

6 

65 

SFC 

18:11:05.69 

-19:37:05.40 

— 

300 

180 

-30 

5.2 

<0.1 

6,2 

72 

HII 

18:11:32.30 

-19:30:42.10 

300 

— 

— 

— 

5.2 

<0.1 

8,6 

73 

SFC 

18:11:46.63 

-18:17:40.40 

— 

240 

120 

90 

3.8 

— 

1 

74 

HII 

18:11:43.26 

-18:16:54.00 

— 

480 

300 

-10 

3.6 

— 

1 

75 

SFC 

18:11:56.19 

-18:48:14.50 

200 

— 

— 

— 

3.2 

— 

1 

78 

SFC 

18:11:59.63 

-19:07:48.50 

400 

— 

— 

— 

4.6 

— 

1 

79 

SFC 

18:12:00.81 

-19:36:00.90 

350 

— 

— 

— 

3.2 

— 

1 

81 

SFC 

18:12:04.07 

-17:52:43.80 

210 

— 

— 

— 

4.3 

— 

6 

82 

SFC 

18:12:08.63 

-16:42:41.60 

360 

— 

— 

— 

3.4 

— 

6 

83 

HII 

18:12:12.42 

-17:40:52.80 

— 

450 

360 

-63 

2.2 

0.4 

4,2,6 

86 

SFC 

18:12:14.29 

-18:26:52.00 

130 

— 

— 

— 

4.8 

— 

6 

88 

SFC 

18:12:26.51 

-17:32:46.30 

— 

240 

90 

-10 

3.9 

— 

6 

90 

HII 

18:12:32.90 

-18:30:13.00 

80 

— 

— 

— 

4.2 

<0.1 

4,3 

91 

SFC 

18:12:32.62 

-17:29:16.80 

— 

120 

100 

80 

3.6 

— 

1 

94 

SFC 

18:12:41.69 

-18:42:58.80 

120 

— 

— 

— 

5 

— 

4 

98 

SFC 

18:12:56.03 

-19:04:23.40 

— 

700 

250 

-80 

3.9 

— 

6 

100 

HII 

18:13:09.20 

-18:07:49.70 

— 

450 

250 

-40 

2.6 

<0.1 

4,2 
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Table R3. Continued. 


ID 

Type 

RA (J2000) 

DE (J2000) 

Radius (") 

Maj. Ax (") 

Min. Ax (") 

Angle (°) 

Distance^ (kpc) 

Dist. dispersion^ (kpc) 

References 

102 

SFC 

18:12:51.72 

-18:48:13.50 

— 

500 

120 

-50 

5 

— 

-6 

103 

HII 

18:13:12.52 

-18:00:07.20 

150 

— 

— 

— 

2.6 

0.1 

4,3 

104 

SFC 

18:13:12.15 

-16:41:11.40 

— 

210 

100 

30 

3 

— 

6 

105 

SFC 

18:13:09.52 

-18:15:57.60 

— 

370 

240 

-45 

3.5 

— 

1 

106 

SFC 

18:13:30.90 

-17:18:06.70 

165 

— 

— 

— 

2.3 

— 

1 

108 

SFC 

18:13:38 

-18:12:15.20 

— 

180 

120 

35 

3.6 

— 

1 

no 

SFC 

18:13:35.59 

-17:23:47.20 

— 

270 

210 

-10 

4.5 

— 

1 

111 

SFC 

18:13:54.20 

-17:16:28.50 

— 

330 

210 

90 

2.4 

— 

1 

114 

HII 

18:13:52.96 

-18:56:44.90 

— 

450 

300 

100 

3.9^ 

0.3 

4,3,8,6 

116 

HII 

18:14:09.82 

-17:27:24.20 

— 

570 

220 

-45 

4.4 

0.1 

4,3,2 

120 

SFC 

18:14:05.27 

-18:14:49.00 

300 

— 

— 

— 

3.5 

— 

1 

121 

SFC 

18:14:13.00 

-18:20:41.80 

— 

150 

120 

-30 

5.1 

— 

6 

122 

HII 

18:14:16.00 

-17:57:18.60 

— 

780 

550 

-35 

4 

0.6 

4,8,3,2 

123 

SFC 

18:14:16.28 

-17:14:57.00 

— 

210 

120 

-60 

2.8 

— 

1 

134 

HII 

18:14:35.87 

-16:45:31.30 

— 

960 

560 

-60 

4.4 

0.3 

4 

136 

HII 

18:14:30.50 

-17:38:34.00 

— 

750 

400 

-65 

4.4 

0.6 

1,4,2,6 

140 

SFC 

18:14:50.37 

-17:20:31.50 

— 

420 

210 

-75 

4.5 

— 

4 

146 

SFC 

18:14:55.69 

-17:46:58.50 

— 

210 

120 

80 

5.1 

— 

6 

149 

SFC 

18:14:57.28 

-18:27:01.10 

270 

— 

— 

— 

5 

— 

6 

152 

SFC 

18:15:01.51 

-17:42:22.80 

100 

— 

— 

— 

4.6 

— 

6 

154 

SFC 

18:15:11.03 

-17:47:21.20 

— 

360 

130 

95 

3.7 

— 

1 

159 

SFC 

18:15:46.15 

-17:34:20.50 

— 

750 

540 

80 

3.8 

0.4 

1,4,2,6 

160 

HII 

18:15:37.51 

-17:04:23.50 

400 

— 

— 

— 

4.1 

<0.1 

4,3 

162 

SFC 

18:15:38.97 

-18:10:42.50 

760 

— 

— 

— 

3.5 

— 

2,6 

163 

HII 

18:15:40.22 

-17:19:46.70 

260 

— 

— 

— 

2.6 

— 

1 

165 

HII 

18:15:45.18 

-16:38:58.20 

240 

— 

— 

— 

5.4 

— 

3 

170 

SFC 

18:16:00.58 

-16:04:53.20 

440 

— 

— 

— 

2.8 

— 

4 

175 

SFC 

18:16:12.68 

-16:43:53.70 

— 

240 

150 

-10 

4.3 

— 

6 

177 

SFC 

18:16:20.00 

-17:05:00.90 

— 

120 

70 

80 

3.7 

— 

1 

178 

HII 

18:16:31.62 

-16:50:59.60 

— 

600 

240 

240 

3.7 

<0.1 

4,2 

179 

HII 

18:16:04.39 

-19:42:05.20 

460 

— 

— 

— 

1.6 

— 

8 

183 

SFC 

18:16:45.65 

-17:02:26.90 

300 

— 

— 

— 

3.7 

— 

1 

185 

HII 

18:16:57.86 

-16:15:33.90 

— 

600 

360 

-60 

2.4 

0.3 

4,3,2,8 

186 

HII 

18:16:51.62 

-18:41:28.30 

300 

— 

— 

— 

4 

— 

4 

187 

HII 

18:16:58.96 

-16:42:08.70 

— 

380 

240 

30 

3.6 

0.2 

4,2,6 

190 

HII 

18:17:11.87 

-16:27:00.60 

— 

630 

300 

-25 

3.T 

0.1 

4,3,2 

192 

SFC 

18:16:18.76 

-14:18:57.50 

1100 

— 

— 

— 

2.2 

0.2 

4 

195 

SFC 

18:17:27.20 

-17:08:58.50 

— 

480 

220 

-45 

2.7 

0.3 

1 

196 

SFC 

18:17:21.85 

-17:00:48.80 

— 

330 

220 

-60 

2.3 

— 

6 

200 

SFC 

18:17:36.02 

-15:56:27.60 

— 

360 

200 

-30 

2.7 

0.1 

4,7 

202 

SFC 

18:17:31.86 

-16:16:59.20 

— 

330 

140 

-80 

3.7 

— 

1 

207 

SFC 

18:17:35.93 

-15:47:45.70 

240 

— 

— 

— 

3.1 

0.2 

1,7 

209 

SFC 

18:17:52.47 

-16:29:02.00 

— 

370 

210 

-45 

4.6 

— 

1 

210 

SFC 

18:17:55.42 

-17:13:05.50 

250 

— 

— 

— 

2.9 

— 

6 

214 

HII 

18:17:45.73 

-16:01:22.20 

— 

400 

200 

-30 

2.8 

<0.1 

7,6 
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Table R3. Continued. 


ID 

Type 

RA (J2000) 

DE (J2000) 

Radius (") 

Maj. Ax (") 

Min. Ax (") 

Angle (°) 

Distance^ (kpc) 

Dist. dispersion^ (kpc) 

References 

217 

SFC 

18:17:59.25 

-16:15:15.20 

— 

480 

270 

100 

3.7 

0.2 

1,2,4 

218 

HII 

18:18:09.85 

-16:52:16.80 

— 

750 

500 

-70 

2.5 

0.4 

1,2,4 

219 

SFC 

18:18:15.15 

-16:04:42.50 

— 

360 

200 

-10 

3.9 

— 

7a 

220 

HII 

18:18:16.55 

-15:59:13.50 

— 

180 

135 

80 

4 

— 

4,3 

227 

SFC 

18:18:46.76 

-16:22:29.30 

— 

170 

150 

80 

4.6 

— 

1 

228 

HII 

18:18:47.08 

-13:32:56.90 

— 

2000 

1200 

-60 

2 

0.4 

4,3,8,2,6 

230 

HII 

18:18:54.20 

-16:48:55.70 

— 

260 

160 

90 

2.7 

0.3 

4,2,8 

233 

SFC 

18:19:06.85 

-16:33:09.50 

— 

570 

280 

100 

2.3 

0.2 

1,4,2 

234 

SFC 

18:19:06.97 

-16:11:20.30 

300 

— 

— 

— 

4.7 

<0.1 

4 

236 

SFC 

18:19:11.92 

-16:18:25.60 

— 

420 

150 

— 

2.6 

— 

1 

238 

SFC 

18:19:38.17 

-16:42:43.10 

— 

650 

330 

-45 

2.4 

— 

1 

242 

HII 

18:19:51.20 

-15:55:05.60 

200 

— 

— 

— 

2.3 

— 

4 

246 

SFC 

18:20:12.44 

-13:51:05.10 

— 

450 

360 

70 

2.6 

<0.1 

7,6 

247 

SFC 

18:20:37.42 

-15:37:55.10 

100 

— 

— 

— 

3.6 

— 

1 

248 

HII 

18:20:27.49 

-16:07:40.60 

850 

— 

— 

— 

2.2 

0.2 

4,8,3,2 

251 

SFC 

18:20:39.33 

-14:02:41.30 

450 

— 

— 

— 

2.6 

— 

1 

253 

SFC 

18:20:56.57 

-15:24:41.80 

180 

— 

— 

— 

3.8 

<0.1 

1,4,6 

256 

HII 

18:21:05.67 

-14:17:49.10 

— 

580 

340 

90 

2.5 

— 

7 

259 

HII 

18:21:09.17 

-14:31:46.20 

— 

700 

350 

-40 

4.3^ 

0.5 

4,3,8 

260 

HII 

18:21:10.92 

-15:02:42.00 

— 

300 

150 

80 

3.3 

— 

3 

261 

SFC 

18:21:12.27 

-16:30:08.20 

210 

— 

— 

— 

2.6 

— 

1 

263 

HII 

18:21:04.64 

_14:46:39.40 

— 

600 

450 

-10 

3.9 

<0.1 

4 

265 

SFC 

18:21:44.83 

-14:56:52.90 

— 

380 

250 

-60 

3.9 

— 

6 

270 

HII 

18:21:59.64 

-14:16:11.60 

— 

360 

260 

-40 

3.6^ 

0.2 

4,7 

271 

HII 

18:22:04.97 

-14:08:55.40 

— 

210 

150 

-30 

3.6 

— 

4,3 

272 

HII 

18:22:05.14 

-14:48:42.90 

100 

— 

— 

— 

3.8 

— 

1 

279 

HII 

18:22:22.95 

-14:35:35.30 

— 

190 

120 

-60 

3.8 

— 

1 

287 

SFC 

18:22:41.01 

-14:27:54.20 

330 

— 

— 

— 

2.4 

— 

7 

292 

HII 

18:22:52.78 

-13:59:26.60 

550 

— 

— 

— 

4.8 

— 

6 

297 

SFC 

18:23:18.85 

-13:15:50.60 

240 

— 

— 

— 

2.2 

<0.1 

5,4,7 

298 

SFC 

18:23:26.37 

-13:49:53.30 

240 

— 

— 

— 

4.3 

— 

7 

301 

HII 

18:23:35.90 

-13:09:28.50 

240 

— 

— 

— 

6.6 

— 

5 

302 

SFC 

18:23:37.57 

-13:18:52.90 

210 

— 

— 

— 

4.2 

— 

5 

313 

HII 

18:24:34.34 

-12:52:02.70 

300 

— 

— 

— 

4 

0.3 

4,3,8 

320 

SFC 

18:25:07.20 

-14:35:12.50 

400 

— 

— 

— 

3.4 

— 

4 

322 

SFC 

18:25:09.19 

-12:44:38.70 

— 

300 

120 

100 

3.5 

<0.1 

1,4 

323 

HII 

18:24:58.83 

-12:01:05.50 

650 

— 

— 

— 

4.5 

— 

6 

325 

HII 

18:25:16.95 

-13:14:13.00 

850 

— 

— 

— 

4.6 

0.5 

4,3,8,2,5 

326 

SFC 

18:25:19.61 

-12:53:00.50 

— 

500 

230 

95 

3.9 

0.2 

1,5,7 

327 

SFC 

18:25:22.32 

-13:34:38.30 

320 

— 

— 

— 

4 

0.5 

1,7 

328 

HII 

18:25:35.01 

-12:20:58.10 

— 

660 

360 

65 

4.2 

0.1 

4,2,5 

330 

SFC 

18:25:28.88 

-13:58:26.00 

420 

— 

— 

— 

5.1 

— 

6 

332 

SFC 

18:25:53.36 

-12:06:08.30 

360 

— 

— 

— 

2.5 

<0.1 

1 

336 

HII 

18:26:00.56 

-11:52:26.90 

— 

720 

440 

-20 

1.9 

— 

4,8,3,5 

340 

HII 

18:26:20.45 

-12:40:52.20 

240 

— 

— 

— 

4.6 

0.3 

4,7,6 
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Table R3. Continued. 


ID 

Type 

RA (J2000) 

DE (J2000) 

Radius (") 

Maj. Ax (") 

Min. Ax (") 

Angle (°) 

Distance^ (kpc) 

Dist. dispersion^ (kpc) 

References 

341 

SFC 

18:26:22.16 

-12:57:16.60 

180 

— 

— 

— 

3.6 

— 

1 

342 

HII 

18:26:24.02 

-12:03:09.10 

— 

200 

160 

80 

2.5 

0.2 

8,3 

343 

SFC 

18:26:30.07 

-12:32:42.00 

— 

360 

300 

10 

4.7 

0.3 

7,6 

348 

SFC 

18:26:46.14 

-12:03:25.70 

230 

— 

— 

— 

4.6 

— 

6 

351 

HII 

18:26:44.64 

-12:24:05.30 

— 

390 

330 

200 

4.5 

0.2 

4,8,5 

359 

HII 

18:27:03.75 

-12:43:48.30 

480 

— 

— 

— 

4.5 

0.3 

5,4,2,3 

363 

SFC 

18:27:17.69 

-10:36:28.60 

300 

— 

— 

— 

2.6 

0.1 

1,7 

365 

SFC 

18:27:09.45 

-13:19:49.80 

480 

— 

— 

— 

4.2 

— 

5 

366 

SFC 

18:27:35.55 

-12:19:14.70 

290 

— 

— 

— 

4 

0.1 

1 

367 

SFC 

18:27:45.09 

-12:46:41.10 

320 

— 

— 

— 

4.4 

0.2 

5,4 

370 

SFC 

18:27:42.08 

-11:36:29.60 

— 

300 

180 

-60 

4.4 

— 

1 

372 

HII 

18:27:52.09 

-12:36:17.50 

440 

— 

— 

— 

4.9 

— 

2,5 

373 

SFC 

18:28:06.14 

-11:39:09.90 

— 

100 

70 

-60 

4.2 

— 

1 

376 

SFC 

18:28:20.95 

-11:35:51.60 

— 

180 

no 

-50 

3.5 

<0.1 

1,4 

380 

SFC 

18:28:23.18 

-11:41:19.00 

— 

260 

220 

-30 

4.4 

0.3 

1,4,7 

381 

SFC 

18:28:23.59 

-11:47:37.90 

240 

— 

— 

— 

3.5 

0.2 

4,5 

386 

SFC 

18:28:48.60 

-11:48:40.00 

— 

420 

230 

-60 

4.6 

— 

1 

389 

SFC 

18:29:18.34 

-12:10:38.60 

— 

820 

430 

-10 

4.2 

— 

5 

390 

SFC 

18:29:14.36 

-11:50:21.70 

250 

— 

— 

— 

3.4 

0.2 

4,2,5 

396 

SFC 

18:29:59.97 

-11:00:30.50 

330 

— 

— 

— 

4.2 

0.2 

1 

399 

SFC 

18:30:30.89 

-11:12:11.90 

— 

630 

330 

-30 

4.7 

— 

5 

400 

SFC 

18:30:30.81 

-11:52:50.40 

730 

— 

— 

— 

3.5 

<0.1 

5,7 

406 

HII 

18:18:47.30 

-15:48:58.30 

220 

— 

— 

— 

1.9 

— 

7 

407 

HII 

18:19:59.95 

-15:27:06.90 

560 

— 

— 

— 

3.7 

— 

6 

26a 

SLC 

18:08:49.40 

-19:40:15.00 

80 

— 

— 

— 

3.1 

<0.1 

2,1 

30a 

SLC 

18:09:08.56 

-19:45:56.20 

60 

— 

— 

— 

5.2 

— 

2 

54c 

SLC 

18:10:40.07 

-19:10:40.90 

100 

— 

— 

— 

3.5 

0.2 

1,3,2,4 

65a 

SLC 

18:11:09.60 

-19:36:32.00 

96 

— 

— 

— 

5.2 

<0.1 

2,6 

83c 

SLC 

18:12:19.40 

-17:40:12.20 

70 

— 

— 

— 

2.2 

0.4 

2,4,6 

106a 

SLC 

18:13:35.58 

-17:18:39.60 

— 

60 

35 

-50 

2.3 

— 

1 

116b 

SLC 

18:14:13.22 

-17:25:15.80 

50 

— 

— 

— 

4.4 

0.1 

3,2,4 

122b 

SLC 

18:13:43.80 

-17:56:34.00 

— 

80 

60 

— 

4 

0.6 

3,2,4,8 

185a 

SLC 

18:16:50.35 

-16:21:44.80 

— 

150 

80 

60 

2.4 

0.3 

3,2,4,8 

203a 

SLC 

18:17:32.59 

-17:06:41.60 

120 

— 

— 

— 

2.6 

— 

2 

217a 

SLC 

18:18:02.99 

-16:17:23.10 

80 

— 

— 

— 

3.7 

0.2 

1,2,4 

218b 

SLC 

18:18:03.23 

-16:54:35.50 

— 

95 

60 

-80 

2.5 

0.4 

1,2,4 

233a 

SLC 

18:19:02.22 

-16:39:59.00 

125 

— 

— 

— 

2.3 

0.2 

1,2,4 

233b 

SLC 

18:19:13.55 

-16:34:57.10 

85 

— 

— 

— 

2.3 

0.2 

1,2,4 

233e 

SLC 

18:19:13.36 

-16:35:04.20 

100 

— 

— 

— 

2.3 

0.2 

1,2,4 

247a 

SLC 

18:20:23.82 

-15:38:31.80 

— 

170 

80 

95 

3.6 

— 

1 

325g 

SLC 

18:25:23.77 

-13:19:04.50 

— 

125 

80 

-10 

4.6 

0.5 

2,3,4,5,8 

325i 

SLC 

18:25:22.56 

-13:17:13.40 

40 

— 

— 

— 

4.6 

0.5 

2,3,4,5,8 

328c 

SLC 

18:25:22.76 

-12:23:26.90 

100 

— 

— 

— 

4.2 

0.1 

2,4,5 

379a 

SLC 

18:28:18.30 

-12:06:26.00 

170 

— 

— 

— 

4.4 

— 

2 

390b 

SLC 

18:29:26.67 

-11:50:45.60 

50 

— 

— 

— 

3.4 

0.2 

2,4,5 
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Table R3. Continued. 


ID 

Type 

RA (J2000) 

DE (J2000) 

Radius (") 

Maj. Ax (") 

Min. Ax (") 

Angle (°) 

Distance^ (kpc) 

Dist. dispersion^ (kpc) 

References 

159a 

SLC 

18:16:06.90 

-17:36:27.10 

65 

— 

— 

— 

3.8 

0.4 

1,2,4,8 

162a 

SLC 

18:15:40.50 

-18:13:08.00 

75 

— 

— 

— 

3.5 

<0.1 

2,6 

160a 

SLC 

18:15:25.30 

-17:05:43.70 

90 

— 

— 

— 

4.1 

<0.1 

3,2,4 

187f 

SLC 

18:16:59.68 

-16:39:26.80 

— 

90 

70 

-40 

4 

0.2 

2,4,6 

212a 

SLC 

18:17:50.60 

-15:53:34.00 

50 

— 

— 

— 

2.7 

— 

2 

212e 

SLC 

18:17:52.84 

-15:55:16.90 

50 

— 

— 

— 

2.7 

— 

2 

253a 

SLC 

18:20:55.50 

-15:24:05.00 

50 

— 

— 

— 

3.8 

<0.1 

1,4,6 

329a 

SLC 

18:25:32.88 

-13:01:51.70 

120 

— 

— 

— 

4.5 

— 

2 

372a 

SLC 

18:27:43.01 

-12:35:42.70 

50 

— 

— 

— 

4.9 

<0.1 

2,5 

397a 

SLC 

18:30:02.90 

-12:15:51.00 

120 

— 

— 

— 

3.1 

— 

2 


Re ferences. (1) iHlsworth-Bo wers et aO (l2013h : (2) iTackenberg et al.l J2012h : (3) lUrquhart et al.l J2013h : (4) IWienen et al.l J2012h : (5) iRoman-Duval et al.l J2009h : (6) iMar shall et al.l (l2009h : 
(7^ ISimon et ^<2006h : rS^ IWalsh et al.Nl997h . 

Notes. These regions have not previous solution for the KDA, but they show shadows against the NIR background radiation, so we located them at the near distance solution. The distance 
shown in this table is the result of averaging all the distance estimates available for each source. The dispersion of all the distance estimates available for each region, when more than one distance 
estimates are available. 
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